Noise effects and tomography of remote entangled spins in quantum dots 
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We investigate decoherent effects on, and propose a tomographic procedure to measure the entangled state 
of two quantum dots in micro cavities, established via the exchange and subsequent measurement of a photon. 
Making the realistic assumption of noise dominated by pure-dephasing processes, we find that the photon must 
be exchanged and measured on timescales shorter than the quantum dots' characteristic dephasing time for 
appreciable levels of entanglement to be achieved. The tomographic scheme is able to reconstruct the full 
density matrix of the quantum dots, and requires only single spin rotations and the injection of an additional 
photon. Remarkably, we find that the additional photon need not be exchanged and measured on a timescale 
shorter than the dephasing time for accurate tomography, and also allows many to be used in order to increase 
the measurement signal. 



I. INTRODUCTION 

Entanglement is a resource in the field of quantum informa- 
tion and plays a particularly important role in quantum com- 
munication and quantum networks!^! As such, there has been 
a huge amount of research into schemes which are able to ef- 
fectively generate entanglement between spatially separated 
systems. Examples include using a direct coupling between 
the two systems; 3 * 5 through projective measurements,^ and 
even via correlated dissipative processesP^l 

One such implementation is the generation of remote spin- 
spin entanglement shared between semiconductor quantum 
dots (QDs) in microcavities via the exchange and subsequent 
measurement of a photonP^El The spin dependent optical 
transitions of a charged QD can cause linearly polarised pho- 
tons to undergo giant Faraday rotations, which can be utilised 
to construct QD-photon entanglement. Allowing a photon en- 
tangled with one such QD to interact with another QD, and 
subsequent measurement of the photon can cause the two QDs 
to become entangled. This scheme is promising since the sub- 
jects of the entanglement - the QDs - constitute good candi- 
date s fo r the storage and manipulation of quantum informa- 
tion! 5 -!^ The photon which transmits the entanglement, on the 
other-hand, is an ideal candidate to transmit quantum infor- 
mation owing to its intrinsically long coherence timelSKD 

In the ideal setting, maximally entangled states of the two 
QDs can be achieved. One source or errors which may cause 
a deviation from this ideal scenario is the unavoidable cou- 
pling of the QDs to their solid-state environments. In typical 
In(Ga)As QDs, the coupling of an excess electron to the nu- 
clear spins in the QD semiconductor via the hyper-fine inter- 
action is thought to dominate over other sources of decoher- 
ence, such as coupling to bulk phonons via spin-orbit inter- 
actions P^KD This can cause a loss of coherence of the elec- 
tron, which will necessarily impair the entangling procedure 
described above. Thus, it is important to establish what lim- 
itations noise is likely to put on the entanglement which can 
be achieved, and to determine how any possible entanglement 
may be experimentally verified. 

In this work we investigate the effects of decoherence on the 
level and type of entanglement which can be obtained. In ad- 



dition, we propose a tomographic scheme which would allow 
for the complete reconstruction of the QD-QD density matrix 
after the entangling procedure has been performed. Our to- 
mographic scheme relies only on single spin rotations of the 
QDs, and the injection of an additional photon. For the pure- 
dephasing noise considered, we find that the tomography is 
unaffected by the time taken for the second photon to be ex- 
changed and measured. As such, the tomographic procedure 
reveals the true state of the QDs after entanglement has been 
established. In addition, we show that it is possible to use 
many photons for the tomography in order to boost the mea- 
surement signal. 



n. SUMMARY OF ENTANGLEMENT SCHEME 

The entanglement scheme first proposed in Ref. [10] relies 
on the spin dependent optical transition from the ground state 
of a singly charged QD to an excited trion state (X~). The 
trion state consists of a pair of electrons each having angular 
momentum projection ±h/2, bound to a heav y hole which has 
angular momentum projection ±2>h/2W&^ There are there- 
fore two excited states, |ti>lT) and -U-), having angular 
momentum projection +3H/2 and — 3h/2 respectively. Sim- 
ple application of the Pauli exclusion principle and conser- 
vation of angular momentum, reveals that an incident photon 
having left-circular polarisation, \L), and angular momentum 
+h, can excite the |^) — > |t-UlT) transition, but does noth- 
ing to Similarly, a right circularly polarised photon \R) 
can excite the \l) — > if) transition, but otherwise does 
nothing. Placing such a QD inside a cavity causes left and 
right circularly polarised photons to acquire different phase 
shifts, depending on the spin orientation of the electron in the 

Following Ref. [10|, we ignore the side leakage of the cav- 
ity k s (or k s < k), where k/2 is the cavity field decay 
rate into the input/output modes, and work in a regime with 
|aj — oj c | <C g and g > (k, 7), where lu and ui c are the fre- 
quencies of external field and cavity mode, g is the coupling 
strength between the trion X~ and the cavity mode and 7/2 
is the X~ dipole decay rate. In doing so, a linearly polarised 
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photon undergoes a giant Faraday rotation with near unity re- 
flectance. The whole process can be expressed as the unitary 

Uiip) = expM\L)(L\ ® ItXtli + IWI ® \m\i)\> (D 

where <p is the difference in phase shifts experienced by left 
and right circular polarisations, and the index indicates inter- 
action of the photon with the ith QD. By adjusting lu — uj c ~ 
±re/2, giant Faraday rotations of tp = ±^ can be achieved!^ 
From a quantum information perspective Eq. ([T]) constitutes 
an entangling gate between the QD and the incident photon. 
By allowing a photon to interact with one such QD in a cavity 
followed by another, an entangled state of the photon and both 
QDs is produced. To see this, we take a vertically polarised 
initial state of the photon, \<f> ph ) = \V) = (1/V2)(\R) + \L)), 
while that of QDs as |<^) = aj|t)i + /3i|4)i for i = 1,2. Thus, 
the initial state of the whole system is \<ptot) — l^ph) (8 ® 
\<j>2), and following the interaction of the photon with both 
QD-cavity systems we have \<p f ) = J7 2 (f )^i(f )|0tot) giv- 
ing 

|^) = |ff>[a 1 a 2 |t)i|t>2-/3i/32|l)i|l) 2 ] 

+i \V)[aM)i\i)2 + "2/3i|l)i|t>2], (2) 

where \H) = (l/v2)(|-R) — \L)) is a horizontally polarised 
photon. Thus, when measuring the photon in the linearly po- 
larised basis, {| V) , \H)} the state of the QDs collapses to the 
(normalised) states 

|$> = v^(aia 2 |t)i|t>2- Pifallhllh) for \H) or (3) 
|*) = V / 2(« 1 /3 2 |t)i|^) 2 + a 2 /3 1 |4>i|t) 2 ) for \V) . (4) 

Upon setting ai j2 and /3i i2 to l/v2, we find |<£>) — > \<j>~) = 

(l/>/2)(|tt) - iw» while |*> = (l/V2)(\U) + 

lit))' which are maximally entangled Bell states. We call 
Eq. ([3]) the outcome and Eq. Q the 'W outcome, alluding 
to symmetry present in the corresponding state of the QDs in 
this idealised case. 



in. NOISE EFFECTS ON THE ENTANGLEMENT 
GENERATION 

Matters relating to the imperfect implementation of the op- 
eration described by Eq. (JTJ have been discussed in Ref. iflOl . 
In this section, our aim is to investigate the non-unitary dy- 
namical evolution of the QDs during the entangling proce- 
dure. We therefore assume that the QD-photon interaction 
happens instantaneously, and consider three time intervals 
during which the QDs undergo decoherent processes: t\ la- 
bels the time from the release of the photon until it reaches the 
first QD; t 2 is the time taken for the photon to travel from the 
first QD to the second, and t% the time from the second interac- 
tion until the photon is measured. During each of these times 
the combined photon-QDs system is assumed to undergo non- 
unitary evolution generated by a master equation of Lindblad 



forrrPl(we set h = 1): 

p = -i[H, p}+^2 (l w l\ - ±{p, Lju^j , (5) 

where p is the density matrix of the combined photon-QDs 
system, H is the free Hamiltonian, and the set {L{\ are Lind- 
blad operators describing the decoherent processes we wish 
to consider. We assume the QDs to be under the influence 
of a magnetic field along the z-axis such that the Hamilto- 
nian for the QDs system reads H = |(<rf + erf), where 
°f = ItXtli - and 6 is tne ne l d strength. We as- 

sume the initial state of the entire system is the separable pure 
state \<ptot) given above, with the coefficients «i. 2 and /3i )2 
set to 1/V2- Following both photon-QD interactions and the 
three periods of time for which the system evolves according 
to Eq. ([5]), the final state of the entire system is 

Pi = c £t3 [U 2 (c ct > [U^ l^totX^totl)^])^] (6) 

where the Liouvillian super-operator is defined to satisfy p = 
Cp in Eq. We then measure out the photon and analyse 
the post-measurement ensemble {pk, Pk}, for fc = $ rep- 
resenting the two measurement outcomes. The probability for 
each outcome is given by pk = Tr[pi(irk ® 1 ® 1)], where 
7r<i> = \H)(H\ and 7i> = |V)(V|, while the post-measurement 
state itself is pk — Tr p h[/0f (|&)(fc| ® 1 <8> l)]/pk, where 1 is 
the identity operator, and Tr p h denotes a trace over the photon 
degrees of freedom only. 

For electron spins in QDs, experiments have measured re- 
laxation times as long as T\ ~ ms! 17 * 23 ! In separate experi- 
ments, spin coherence decay timescales have been measured 
to be T 2 ~ ps, at best! 2 ^ As such, we expect any superposi- 
tion of |t) and \l) - as is required for the QD-entanglement 
scheme - to be affected by pure-dephasing processes on a 
timescale T 2 , well before spin relaxation processes become 
important. We therefore consider the pure-dephasing form of 
Eq. |5]), which is achieved with the pair of Lindblad operators 
Li = yr 2 /2of for i = 1, 2, and where T 2 = 1/T 2 . 

We first remark on two important features. Firstly, the mea- 
surement probabilities are unaffected by the non-unitary evo- 
lution of the QDs; the probability that the post-measurement 
state is p$ or p^ remains the same and equal to 1/2 for all 
times. This is perhaps to be expected, since pure-depasing 
does not affect the diagonal elements of the QD density matrix 
when expressed in the basis {|f) , \i)}, and it is these states 
through which the QDs are coupled to the photon. Put an- 
other way, although there is an asymmetry between the dots 
owing to the fact that the photon interacts with one first, this 
asymmetry cannot be distinguished by pure-dephasing noise 
and the form of QD-photon interaction. Our second remark, 
and another manifestation of these symmetry arguments, is 
that the post-measurement states depend only on the sum of 
the three time intervals, t = t\ + t 2 + t 3 . Thus, for pure- 
dephasing noise considered here, there is no particular waiting 
period which affects the post-measurement states more than 
the others. 

The entanglement fidelity of the two post-measurement 
states is defined as — J (fc| pk \k) , for k — $, ^, where 
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FIG. 1: Entanglement (solid red) and entanglement fidelity of the ^ 
(blue dashed) and $ (blue dot-dashed) post-measurement states of 
the QDs, as a function of the total time t — ti + t 2 + £3 of pure 
dephasing. For the entanglement fidelity of the $ state we have set 
e = 10/Ta. 



= 1^ ) anc * l^)o = IV'i) ^ the states obtained in the 



absence of noise. From Eq. (|6| we find 
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FIG. 2: Summary of method for QD tomography. The first photon 
interacts with the QDs and then is measured, leaving the QDs in 
one of two post-measurement states, which we describe by a set of 
coefficients oiij [see Eq. (10) ]. We then immediately apply single 
spin rotations to the QDs. Next, a second photon interacts with one or 
both of the quantum dots. The measurement statistics of the second 
photon then reveal one of the oiij coefficients, depending on which 
rotations were performed on the QDs. 



while 



\0 



-2t/T 2 ) 



(8) 



and where the cosine factor appearing for F$ is present since 
the magnetic field in the z-direction causes oscillations be- 
tween the entangled states \<j)~) and \(j> + ) = (l/-\/2)(|tt) + 
|44) )■ We quantify the entanglement itself using the entangle- 
ment of formation, 25 given by, 

E{t) = -x(t) log 2 [z(*)] ~ (1 - x(t)) log 2 [l - x(t)] (9) 



defined in terms of the function x(t) 



(1/2)(1 



WT— C(t) 2 ) which itself depends on the concurrence C(t), 
which in the present case takes on the simple form C(t) — 
e _2 */ T2 . We find that both the entanglement and the en- 
tanglement fidelity decrease in an exponential fashion with 
time r* 26 * 27 ! The non-zero value of fidelity approached reflects 
that the large t limit of the QD system is the maximally mixed 
state which has non-zero overlap with the target states \4>~) 
and 

In Fig. ([TJ we plot entanglement (red solid curve) and the 
entanglement fidelities (blue dashed) and Fq, (blue dot- 
dashed) as a function of t = t\ + t 2 + t 3 in units of the pure- 
dephasing timescale T 2 . In order for greater levels of entan- 
glement to be reached, all sequence operation times must be 
kept as small as possible: we require that ti,t2,ts < J2. 



IV. TOMOGRAPHY OF QD STATE 

We now propose a simple method to perform tomography 
of the post-measurement state of the QDs. The general idea is 



to apply only single qubit rotations to the QDs, to then inject 
a second photon into the system, and to determine the state 
of the QDs through measurement of the second photon. The 
great advantage to using only single spin rotations being that 
they are more easily achieved experimentally, and that they 
cannot affect the level of entanglement shared between the 
QDs. The tomographic scheme is depicted in Fig. 

Recent experiments have demonstrated rotations of elec- 
trons in QDs on picosecond timescales. 28 We therefore as- 
sume that the rotations we require are effectively instanta- 
neous compared to our other timescales of interest, namely 
T 2 ~ /is and Ti ~ ms. We begin by writing a completely 
general density matrix for the QDs in its Hilbert-Schmidt de- 
composition: 



PQDs - 

where we define aoo 
the real coefficients a; 



1 3 

i,j=0 



1 and cr 



(10) 



We note that 



T^iPQDsCi ® <Jj) are expecta- 
tion values of measurements made on the QDs, and are a 



complete representation of the state. Eq. ( 10 1 represents the 
post-measurement state of the QDs, i.e. p<§, or p^. Follow- 
ing the measurement, we immediately apply single spin rota- 
tions to the QDs with general unitary transformations of the 
form Ri(0) — exp[|0 • erj, where er^ = (af , of , of), and 
6 = (8 X , y ,9 z ) is a vector of the rotation angles. We write 
the state of the QDs after the rotations have been performed as 
Xqds = Ri(ei)R 2 (e 2 ) PQDs R 2 (e 2 yR 1 (e 1 )K where 0i and 
6 2 describe the rotations performed on the two QDs 

The second photon injected into the system is also verti- 
cally polarised, giving a complete density operator p tot — 
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|V}(V| <E> XQDs before the second photon is injected. Once 
the second photon has been allowed to interact with the QDs, 
we trace out the QD degrees of freedom to obtain the reduced 
density operator of the second photon only. This reduced den- 
sity operator is given by 



Tr, 



QDs 



(11) 

where TrqD S is a trace over the degrees of freedom of the 
QDs, and the times Tj for j — 1,2,3 parameterise the second 
photon's path through the QD system. The (1,2) superscript 
indicates that the photon interacts with both QDs. Regardless 
of the rotation angles, we find the state of the second photon 
has the general form 



1 -£ 
-t 1 



(12) 



where the matrix is written in the circularly polarised 
{ \R) , \L}} basis. The real quantity £ is in general a function of 
the set of coefficients ay, the rotation angles Q\ and 62, and 
the times ta. Interestingly, for pure-dephasing noise, however, 
we find that £ does not depend on any of the time intervals 73 
describing the second photon's path through the QD system. 
Thus, any pure-dephasing noise affecting the QDs after the 
necessary rotations have been applied to them, does not affect 
the accuracy of the tomographic scheme. This can be appre- 
ciated by noting that U\ and U% couple only to the diagonal 
elements of xqd S i which are unaffected by pure-dephasing 
noise. As such, the trace in Eq. (Ill picks out density matrix 
elements which are protected from noise. 

We find that with various combinations of the rotations 
X = exp(— Ifffx) and Y — exp(|f a y ) (and the identity) 
on dots 1 and 2, we obtain a reduced density operator of the 
second photon of the form Eq. ( 12 1 where £ is a single value 



ay with i,j ^ 0. Therefore, measurement of the second pho- 
ton in the usual linearly polarised basis, {\V) , \H)}, would 
then lead to outcome probabilities of the form Pq> = i(l — £) 
and Pj> = 1(1 + £) for the ^ and $ outcomes, respectively. 
Thus, analysing the statistics of this measurement it is possi- 
ble to extract the value of £ from £ = Pj> — Pj, and, thus, 
information about the state of the QDs. In Tab. we tabulate 
which rotations, and on which QD(s), give which value of ay. 

From Tab. <|TJ> we see that applying no rotations to the QDs 
results in the second photon containing information about the 
expectation value a zz = Tr(pQ£, s a z <g> a z ). This can be at- 
tributed to the form of the QD-photon interaction, Eq. ([T]). The 
QD-photon interaction depends on the spin projection of the 
QD along the z-axis. Thus, a photon having interacted with 
two QDs in such a way, will have a state whose form depends 
on the correlation of the spin projections, a zz . Rotations ap- 
plied to the QDs amount to a rotation of the matrix with ele- 
ments ay. The rotations therefore move a new element into 
the position which was previously a zz . It is then this new el- 
ement which is protected from pure-dephasing noise, and this 
element which the photon state then depends on. 

We also note the useful property that the state of the QDs 
following measurement of the photon collapses to one hav- 



TABLE I: Those ay coefficients that can be obtained from the re- 
duced density matrix of the second photon by applying rotation(s) to 
the QDs. 



QDi 


QD 2 




1 


1 - 


-> Ozz 


1 


X - 


&zy 


1 


Y - 


-> Ctzx 


X 


1 - 


a yz 


Y 


1 - 


-> a xz 


X 


Y - 




Y 


X - 


' O^xy 


X 


X - 


— > Olyy 


Y 


Y - 





ing a zz = 1 or a zz = — 1, depending on the two measure- 
ment outcomes <I> and ^ respectively. Therefore, if an addi- 
tional photon is injected (after the second) without perform- 
ing any rotations of the QDs, since Pj> = 4(1 — a zz ) and 
P$ = |(1 + a zz ), this additional photon will be measured 
having the same polarisation as the second with unit proba- 
bility. This can be done continuously to build up a string of 
photons all of which will be measured having the same po- 
larisation, with the probability that this polarisation is H or 
V corresponding precisely to the value ay determined by the 
rotations performed after the first entangling photon. 

In order to obtain the remaining a^o and aoj coefficients, 
we allow the second photon sent through to interact with only 
one of the two QDs. In place of Eq. 



1 1 1 we then have 



(2) 

with a similar expression for p\ . In this case we obtain a 
reduced density matrix for the second photon of the form 



Ptot 



<ul] 



(13) 




(14) 



(2) 

with an identical expression for p\ ' , and where £ is a real 
constant. For the same reasons as before, we find that for pure- 
dephasing noise ( does not depend on the times Tj . Applying 
again tt/2 rotations with X and Y, we find that ( is equal 
to one of the remaining coefficients [see Tab. ([TTJ>] . To extract 
the value of £ from the photon, we now need to measure the 
photon in the diagonal | + 45°) = (\R) + i\L))/y/2 and | - 
45°) = (\R) — i\L))/y/2 basis. The outcome probabilities 
would then be P +45 o = \ (1 - Q and P_ 45 o = \ (1 + £)> and C 
can be found from the quantity — P+ 450 . As before, the 

post measurement state of the QDs is such that any additional 
photons injected into the system will be measured having the 
same polarisation with unit probability, and can therefore be 
used to strengthen the measurement signal. 

Thus, we see that with only single spin rotations of the QDs, 
we are able to reconstruct the complete post-measurement 
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TABLE II: Coefficients and aw that can be obtained from the 
reduced density matrix of the second photon by applying a rotation 
to one QD before the second photon interacts only with one of the 
QDs. 



QDi 


QD 2 


c 


1 


(no interaction) — 




X 


(no interaction) — 


-» OLyQ 


Y 


(no interaction) — 


-> QiO 


(no interaction) 


1 - 


-> «0z 


(no interaction) 


X 


-> a 0y 


(no interaction) 


Y - 





state. Moreover, owing to the form of the QD-photon inter- 
action, pure-dephasing noise affecting the QDs after the rota- 
tions have been preformed has no bearing on the accuracy of 
the tomographic procedure. Thus, while it is important for the 
initial photon's path through the QD system to be achieved 
on a timescale shorter than T2 (in order to generate an ap- 
preciable amount of entanglement), the amount of entangle- 
ment that is measured by the second photon reflects the true 
amount that was present immediately after the rotations are 
performed. Additionally, any photons injected after measure- 
ment of the second will be measured having the same polar- 
isation as the second with unit probability, and can therefore 
be used to strengthen the signal. 

It should be noted that the arguments above hold only for 
pure-dephasing noise. Spin relaxation does affect the diago- 
nal elements of the density operator of the QDs. As such, if 
the parameters describing the second photon's path through 
the QD system, Tj, were of the order of the spin relaxation 
timescale, Ti, we would find that £ and £ do depend on Tj. 
Thus, for accurate tomography to be achieved, it is still nec- 
essary that the second photon traverses the QD system on a 
timescale <C T\ . 



V. SPIN-RELAXATION 

Though it seems that with current technologies the dephas- 
ing timescale of electron spins in QDs is likely to be far shorter 
than the spin-relaxation timescale, it is nevertheless interest- 
ing to briefly investigate what differences may occur if the 
coherence of the QDs were limited by relaxation processes. 
In order to do so, we now consider the set of Lindblad opera- 
tors given by Li = \fT\a~ for i = 1,2, where = |4)(t|j 
and T^ 1 = T\ is the spin-relaxation time. 

Our current protocol consists of the injection and measure- 
ment of a first photon to establish entanglement, the appli- 
cation of rotations to the QDs, followed by the injection and 
measurement of a string of photons whose polarisations re- 
veal information about the QD state established. For noise 
originating from spin relaxation processes, the first stage of 
this protocol is largely unaffected. For appreciable levels of 
entanglement to be achieved, the first photon's path through 



the QD system, characterised by the times tj, must be short, 
but now compared to the spin-relaxation time T\ . We do note, 
however, that now the measurement outcome probabilities do 
change, and we have 

= ^e~ tiri (l + e~ t2ri - e - (tl+t2)ri ) , (15) 

and P$ = 1 — Py. Thus, in this case we see that the three time 
intervals play a decreasingly important role in the evolution 
of the probabilities. We also see that as t\ and t<i increase 
the probability of obtaining the $ outcome tends to one while 
that of the ^ outcome tends to unity. This can be understood 
by noting that spin relaxation transfers population to the state 
|44-) which has overlap only with |$). 

For spin-relaxation we find that the tomographic proce- 
dure described in the previous section still works, though 
does not share some of the helpful features previously de- 
scribed. Namely, while a reduced photon state of the form of 
Eq. ( fT2} is found, for the combinations of rotations described 
in Tab. the value £ depends on a combination of ay val- 
ues, as well as the times describing the second photon's path 
through the system. We note, however, that tomography ought 
still to be possible, but will require specific knowledge of Ti, 
and the ability to vary the time taken for the second photon to 
pass through the system. 

We lastly note an intriguing feature of spin-relaxation noise, 
which is that it allows for the level of entanglement to be 
boosted by additional photons which pass through the system 
and measured. Numerical simulations suggest that if photons 
can be sent through the system and measured at a sufficient 
rate, sequences of measurement outcomes all being \& will 
correspond to a state of the QDs with entanglement main- 
tained at a particular level. We believe future research in this 
direction to be worthwhile, though beyond the scope of this 
work. 



VI. SUMMARY 

In Ref. ifTUl a scheme to entangle two spatially separated 
QDs in micro cavities was proposed. The scheme relies on 
a QD-state dependent rotation of the polarisation of a photon 
which is exchanged. In this work we have investigated the 
effects of non-unitary dynamical evolution of the QDs, caused 
by the coupling to their solid-state environment. We found 
that pure-dephasing of the QDs necessarily decreases the level 
of entanglement that can be attained. 

We then proposed a method to perform tomography of 
the state of the QDs, which relies on only single spin lo- 
cal rotations of the QDs, and the injection of additional pho- 
tons. Interestingly, while the level of entanglement attained 
is sensitive to pure-dephasing noise, the accuracy of the to- 
mographic procedure is not. Thus, the time in which to per- 
form the tomography is limited only by the spin relaxation 
timescale, which is typically orders of magnitude greater than 
the pure-dephasing time. Lastly, we found that within the 
pure-dephasing time, many photons can be injected into the 
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system in order to boost the tomographic measurement sig- 
nal. 
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